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Abstract 

A general formula for the entropy of stationary black holes in Lovelock gravity 
theories is obtained by integrating the first law of black hole mechanics, which is 
derived by Hamiltonian methods. The entropy is not simply one quarter of the 
surface area of the horizon, but also includes a sum of intrinsic curvature invariants 
integrated over a cross section of the horizon. 



^ j acobson@umdhep 

^ rcm@hep .physics . mcgill. ca 



The study of black hole thermodynamics is motivated primarily by the hope of learning 
something about the nature of quantum gravity. Probing the limits of validity of the four 
laws of black hole thermodynamics [|I|] may provide one source of insight in this context. For 
instance, one would like to know to what extent these laws are valid when the back-reaction 
of quantum field energy is taken into account. Back-reaction leads to the consideration 
of higher curvature interactions, which arise from quantum field renormalization[0. This 
observation motivates the investigation of black hole thermodynamics in higher-curvature 
theories of gravity. Another motivation comes from string theory, where the low energy 
effective field theory of gravity contains higher-curvature terms ||^. 

In this letter, we focus on a special class of higher-curvature theories, called Lovelock 
gravity^. These theories are the most general second-order gravity theories in higher di- 
mensional spacetimes. Higher dimensional spacetimes are of interest in several candidate 
frameworks for unifying gravity with other interactions. Moreover in higher dimensions 
with Lovelock theory, one can explore the effect of higher curvature terms in black hole 
thermodynamics without having to deal with complications that arise in true higher deriva- 
tive theories. The Lagrangian density for Lovelock gravity in a spacetime of dimension D 
may be written C = J^ml^'^m^m where 

r ( „\ ^ I ~ raibl - Qmbm p C\d\ T3 Cmdm / 1 ^ 

and Cq = \/—g^, and in the sum, [D/2] indicates the integer part of D/2. The 6 symbol 
is a totally antisymmetric product of 2m Kronecker deltas, normalized to take values 
and ±1. Note Ci = ^^—gR by itself yields the Einstein Lagrangian. In general. Cm is the 
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Euler class for a 2m-dimensional manifold. Because of the antisymmetrization, no deriva- 
tive appears at higher than second order in the equations of motion. Static, spherically 
symmetric black hole solutions have been found for Lovelock gravity]^, and given these 
solutions, one can examine to what extent the laws of black hole thermodynamics apply |Q. 

Using the Euclidean approach^, it is clear many of the essential features survive. Given 
a specific stationary solution, a Euclidean section is obtained by analytic continuation. 
Regularity of this section requires identifying the Killing time coordinate with a period (3 = 
271 /k. This section is regarded as a background in a periodic Euclidean path integral, which 
is interpreted as giving a thermodynamic partition function Z with inverse temperature, f3. 
In a semi-classical approximation, the Euclidean action I{f3, J, ■ ■ ■), regarded as a function 
of (3 and the extensive parameters (such as the angular momentum, J), is identified with 
— In Z. One defines the internal energy and entropy as 

dl dl 

"-TP 

which then automatically satisfy the first law 

This result reproduces the first law of black hole mechanics iiU = M, the black hole mass, 
and {Tdl/dJ) = fl, the angular velocity of the horizon. The latter relations appear to 
hold in the present context {i.e., higher curvatures and/or higher dimensions), and it may 
be possible to prove them by extending the methods of Ref. [|l^. One might expect that 



S coincides with one quarter the surface area of the horizon |Tl| as in Einstein gravity, but 



this identity fails in Lovelock gravity 0, and other higher curvature theories . 

As described above, the Euclidean approach is applied to investigate known stationary 
solutions. While this yields S" as a function of the solution's parameters, it provides no 
guidance as to any geometric significance of the black hole entropy. It also fails to provide a 
definition for 5* that generalizes to an arbitrary, time-dependent horizon, without which it is 
difficult to see how to even address the question of whether a classical second law is obeyed 
in "nonequilibrium" processes. However, it is also possible to work with the Euclidean 
approach without specific solutions in mind|TE|, [10]. (Indeed, this has recently been done 



for static metrics in Lovelock gravity with only ci and C2 non- vanishing When applied 
in this more general fashion, it should be possible to avoid the abovementioned drawbacks 
of the Euclidean approach. 

Our alternate approach is to find the variation 6S of the entropy by deriving a general 
form of the first law of black hole mechanics following the Hamiltonian method introduced 



by Sudarsky and Wald||15[. A priori, there is no guarantee that the resulting 6S is the 
variation of a functional of the intrinsic geometry of the horizon, however, one finds that 
it is for Lovelock gravity. In fact, the entropy is given by an integral over a cross section 
of the horizon, of an integrand which is similar in form to the Lagrangian itself. 

The method of Ref. applies to Killing horizons. A Killing horizon is a null hyper- 
surface whose null generators are tangent to a Killing field. In four-dimensional Einstein 
gravity. Hawking proved that a stationary black hole event horizon must be a Killing 



horizon||l6l. This proof can not obviously be generalized to Lovelock gravity, but the result 
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is certainly true for all the known solutions (which, however, are all static). 

If x"" is the Killing field which is null on the horizon, the surface gravity k can be 
defined by V aix'^Xb) = — 2KXa, or equivalently, x^^bX"" = i^X""- If oiie assumes (as we will) 
that the null generators of the horizon can be extended to the past as complete geodesies, 
three important properties follow: (i) k is constant over the entire horizon [T^ . (ii) There 
exists a (i5-2)- dimensional spacelike cross-section B of the horizon on which x"" vanishes 
r?| . i? is called the bifurcation surface, (iii) The extrinsic curvature of B vanishes |I8| . 



A further property of Killing horizons that is important for our purposes is that any two 
spacelike slices of the horizon are isometric. 

Let be the Killing field that is asymptotically a pure time translation. Then by 
suitably rescaling x", one is able to reduce the Killing field x"" — ^'^ to a rotation. In 
higher dimensions, there are [{D — l)/2] commuting Killing fields generating rotations 



in totally orthogonal planes p9|] , and so one has 

x'^ = r + ^^^"'C) 



Next, we briefly present the Hamiltonian description of Lovelock gravity EO], which will 



be necessary in the following derivation of the first law. As usual after splitting space and 
time, the dynamical variables on a spatial surface E are the spatial metric hab, and its 
conjugate momentum 7r"''[^. In the present case, tt°'^ is a complicated function of both the 
extrinsic and intrinsic curvatures of The normal and spatial components of the time 

flow vector field t'*, the lapse N and shift A^", are arbitrary Lagrange multipliers in the 
Hamiltonian. For asymptotically fiat space, the Hamiltonian has two parts H = Hy + Hs- 



The volume term Hy is a combination of constraints 1-L± and Tia, 

Hv= I d^-^x {NH^ + N^Ha) , 
Jt. 

and hence vanishes when evaluated for solutions of the field equations. Tia is the generator 
of spatial diffeomorphisms in S, and so, as in Einstein gravity, Tia = —2Dh{TTj'), where Da 
is the covariant derivative compatible with hab- 'H± generates normal deformations of S, 
and one finds 7{± = Scm'^^i"^ with 

-lyC™-) ^aibi---am.bm JD cidi p Cmdm 

'^1. ~ 2m "cidi---Cmdm -'^aibi ' ' ' -'^ambm 

where Rab'^'^ is the full D-dimensional curvature tensor. Since Rab'^'^ is projected into S, it 
can be replaced by Rab'^'^ + '^K^a'^Kf,]'^. Here Rab'^'^ is the curvature of hab, and Kat is the 
extrinsic curvature of S, which is regarded as a function of hab and 7r''^||2T|. The 6 symbol 
is the antisymmetric product of projected Kronecker deltas, 6^ = 6^ + n°'nc where n° is the 
unit normal to S. 



The Hamiltonian also has an asymptotic boundary term, just as in Einstein gravityp2|, 
Hs = §ao d^~'^x {NS± + N°-Sa)- Hs is added to cancel surface terms which arise by integrat- 
ing by parts in Hy when generating the Hamiltonian flow equations. The precise value of 
Hs depends on the choice of = Nn°' + N°'. If asymptotically C,"', the time translation 
Killing field, then Hs = M, the total mass of the solution. If — (j)'^^^, a rotational Killing 
field, then Hs = —J(a), the associated angular momentum. 



The derivation of the first law comes from a judicious evaluation of the Hamiltonian 115 
Let hab and vr"^ describe a stationary black hole solution, and choose the "time" flow field 



to be the Killing field that is null on the horizon, = = + ^'"Vfa)- Evaluate H on 
a space-like slice E that extends from asymptotic infinity to the horizon, intersecting the 
horizon at the bifurcation surface, B. With these choices one has H — M — O*^"' J(„). Let 
{Shab, Stt"''') be a perturbation of the initial solution to any nearby solution — not necessarily 
stationary. With and E fixed, one has 

5H^5M- Q(")5J(«) 

On the other hand, one has Hamilton's equations 

The volume term vanishes because the flow is along a KiUing field, C^hab = = >C^7r"^. 
Integration by parts is needed to produce these volume terms. The surface terms at 
asymptotic infinity are cancelled by the variation 5Hs, but one is left with surface terms, 
denoted SHb, at the inner boundary of E. 

Since SHb is evaluated on the bifurcation surface B where t"" = vanishes, any 
nonvanishing terms involve derivatives of t"'. Such terms only arise from the metric 
variation of the curvatures in N7i±. Two integrations by parts arise from SRab^'^ — 
-2DyaDHhb]'^ + Rab^HhJ^ , where Sha'^ = hl^^Shab. The first integration yields new volume 
integrals involving DaN . In these integrals, the second integration by parts produces a 
boundary term at B involving DaN. The contribution from Ti^"* is 

^^b'"^ -^if/^'-'x [^hv^Wa,NS^^l'/^Z^^^^^^ (3) 
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where v'^ is the unit normal to the bifurcation surface B, pointing into S, and h is the deter- 
minant of hab, the induced metric on B. Now it is not hard to show that, at B, DaN = nva, 
where k, is the surface gravity. Since Vai^^^^cldl-'-dZ ~ ^dl '- ^dl^ antisymmetrized product 
of Kronecker deltas S'^ = S'^ + n°-nc — v°'Vc projected into B, all remaining tensors are fully 
projected into B. The projected Shb"^ yields 6hb'^ = h'^'^dhbe- Since the extrinsic curvature 
of the bifurcation surface B vanishes, projecting Rab^'^ yields the intrinsic curvature Rab^'^ 
of Kb- 

Thus far variations of the Kab terms in Rab'^'^ were neglected. Regarded as a function of 
tt"* and hab, Kab involves Rab'^'^ and so could produce extra boundary contributions. These 
extra terms always vanish however, since they contain at least one other factor of Kab 
which vanishes when projected into B. Since k is a constant on S, Eq. (^) thus becomes 

Extracting the temperature, k/2ti, the contribution 55'^™'' = {2ti / k)6H'^^'^ to the variation 
in the entropy thus depends only on the intrinsic geometry of B. The next step is to 
recognize that is in fact the variation of an intrinsic geometric quantity on i?, S**^"^-* = 

Aumcm §3^^'"^ ^m-i{h) with Cm as defined in Eq.(|l|). 

Collecting all of the contributions to the entropy, our final result is the first law of black 
hole mechanics in Lovelock gravity, 5M — f2*°'(5J(a) = {k,/27t)6S where 

S=j2 47rmcmf d''-^ Cm-i{h) (4) 

m=l 

Our derivation of the first law relied heavily on the fact that the variation of the entropy 
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is evaluated on the bifurcation surface B. Nevertheless, the integrated expression (^) for 
the entropy of a stationary black hole can be evaluated on any spacelike slice of the Killing 
horizon, since all such slices are isometric. Note that with ci = I/IGttG, as appropriate 
for the Einstein action, the first term in the sum is simply A/ {AG), where A is the surface 
area of the horizon. It is remarkable that the entropy @ is identical in form to the original 
action, evaluated in a (Riemannian) space of dimension D — 2, with Cm-i replaced by 

For even dimensions, we have included a non-trivial integration constant in Eq. (^). 
This constant is the contribution of £d/2-i, which yields the Euler constant of the cross 
section of the horizon. In four dimensions, this constant is fixed for all stationary black 
holes since the horizon must have the topology of S^xi?. |16| In D>4, the horizon topology 



is not unique, so this topological entropy may have more significance. With this choice for 
integration constant one easily confirms that for the known black hole solutions 0, Eq. 
produces the same result as Eq. (||) in the Euclidean approach (provided, of course, that 
cd/2 is identical for both approaches). 

The present derivation is easily extended to include Maxwell or Yang-Mills fields [p!5|. 
In the case of charged black holes, the first law becomes 5M — — f2'"'5J(„) = {n/27!-)6S, 
where $ is the potential difference between the horizon and infinity, and S is still given 
by Eq. @. Again, Eq.'s (H) and (j^) yield the identical results when applied to the known 
charged black hole solutions 0. 

Our results extend the framework of black hole thermodynamics in a natural way to 
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Lovelock gravity. The zeroth law {i.e., the constancy of k) holds modulo the assumption 
that the horizon generators are geodesically complete (or, if the black hole forms from 
collapse, that the horizon generators of the stationary solution which it approaches have 
this property). The first law defines an entropy that is localized in the intrinsic geometry of 
the horizon and is valid for all stationary black holes, at least provided the event horizon is 
a Killing horizon. For example, Eq. applies for rotating black hole solutions, which are 
as yet unknown. If all the coefficients Cm are of the same order in units of some common 
length scale, then the area term (m = 1) dominates the entropy for black holes that are 
much larger than that scale. The entropy of small black holes can be negative, but for 
fixed couplings Cm it is bounded below for the known solutions|P]. Such negative values 
are thermodynamically benign, since only changes of the entropy occur in the first law. 
The second law for quasi-stationary processes follows from the first law, provided the 



quantity 5M — i7*"'(5J(^) is positive. This is the case for fiuxes of positive energy matter p3 
To control the sign for fiuxes including gravitational energy would require a positive energy 
theorem for Lovelock gravity, which has not been established. p4[ (Negative energy states 
would probably make the theory unstable however, in which case a failure of the second 
law would be of less significance.) Alternatively, the second law in Einstein gravity can be 



derived from Hawking's area theorem^] which shows that, even in non-quasi-stationary 



processes, the area (and therefore the entropy) will never decrease provided the Ricci tensor 
(or, using Einstein's equation, the matter energy tensor) obeys the null energy condition. 
In Lovelock gravity the expression (HI) for the entropy is at least meaningful on a slice of 
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a time-dependent horizon, so one might hope to prove a non-equihbrium second law by 
similarly following the evolution of the horizon. This remains an interesting problem for 
future work. In this connection curiously, note that with 0^/2 > the topological term in 
Eq. (0) could lead to violations of the second law when two black holes coalesce, even in 
four dimensions. 

When quantum fields are included, negative energy can be transferred to the black 



hole, as in Hawking evaporation p5| or "mining" processes. The interesting question is 
then whether the generalized second law (GSL) {S{Sbh + Soutside) > 0) holds. Validity of 
the classical second law for 5*^^ alone would seem to be a prerequisite for validity of the 
GSL. Provided Lovelock gravity turns out to satisfy the classical second law, it seems to 
be on the same footing with Einstein gravity as far as the GSL is concerned. For instance. 



Zurek's argument ||27|| that entropy increases during Hawking evaporation uses only the first 



law, and not the form of the black hole entropy, so it will carry over to Lovelock gravity. 



More general arguments]^, to the extent that they are valid, also seem to carry over. 
However, even if valid, these arguments apply only to quasi-stationary processes, so the 
validity of the GSL remains an important open problem. 

Finally, Lovelock gravity is only a very special case of possible generally covariant 
gravity theories. The fact that the entropy is not one quarter the surface already 



known in many other examples of theories with higher curvature interactions [|T2|. If one 
derives a first law via the method of Sudarsky and Wald as we have done here, the variation 
of the entropy will again be given as an integral over the bifurcation surface of the horizon. 
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However, this expression will not in general depend only on the intrinsic geometry of the 
horizon, and it is not even clear that it will be the variation of some quantity defined locally 



at the horizon (although, from the Euclidean approach, Ref. [14] argues that the entropy 
is always a quantity defined locally at the horizon). In fact, the particular properties of 
Lovelock gravity played a crucial role in our calculations establishing the local and intrinsic 
nature of the entropy. Whichever way it turns out in general higher curvature theories, 
one of the authors of this paper will enjoy a free dinner at a restaurant of his choice. 
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